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Abstract Recently, de Roany & Pacheco [1] performed a Newtonian analysis on the 
evolution of perturbations for a class of relativistic cosmological models with Creation 
of Cold Dark Matter (CCDM) proposed by the present authors [2]. In this note we 
demonstrate that the basic equations adopted in their work do not recover the specific 
(unperturbed) CCDM model. Unlike to what happens in the original CCDM cosmology, 
their basic conclusions refer to a decelerating cosmological model in which there is no 
transition from a decelerating to an accelerating regime as required by SNe type la 
and complementary observations. 
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An alternative cosmological approach providing a transition from an early decel- 
erating to a late time accelerating expanding Universe (as indicated by SN la data 
[2) and a reduction of the so-called cosmic dark sector has been recently discussed 
in the literature [5]|31[S]. The basic idea is a simple one. The gravitationally- induced 
particle creation can lead to an accelerating cold dark matter (CDM) dominated Uni- 
verse without requiring the presence of quintessence scalar fields or a cosmological 
constant. This happens because the irreversible creation of cold dark matter (CCDM) 
can thermodynamically be described by a negative pressure thereby resulting in a pos- 
itive acceleration as predicted by the Einstein field equations [6]. By neglecting the 
radiation and baryonic components, a pure CCDM cosmology can be fully described 
by the following relativistic equations [2] 



A 

, a k 



(1) 



8nGpc = -2 n-^, (2) 

where p is the CDM density, pc is the creation pressure, and an overdot means time 
derivative. In the case of constant specific entropy per particle ("adiabatic" particle 
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creation), the creation pressure for a CDM component is given by [6] 

Pc--§, (3) 

wiiere H = d/a is the Hubble parameter and F is the creation rate of CDM particles. 

Now, following standard lines, the above equations can be rewritten in order to 
obtain the "continuity" and acceleration equations: 

p + 3Hp = pr, (4) 

-{p + 3pc) = 1-T7 • (5) 



The above Eq. ((Sjl shows that creation of CDM particles [F > 0) may provide a 
transition from a decelerating to an accelerating stage. Note also that the relativistic 
description of a CCDM cosmology is fully determined once the F parameter is given. In 
the specific CCDM cosmology proposed by Lima et al. [2] (from now on LJO model), 
the creation rate is defined by: 

;^=3a— , (6) 

where pco = SHq /SttG is the present value of the critical density. Inserting the above 
expression into one finds for the creation pressure, Pc = —apco, which is negative 
and constant (termed -A in the notation of [1]). At the level of background equations, 
it has also been proved that the cosmic history of such a CCDM cosmology is indistin- 
guishable from the standard ylCDM model. As one may check, by combining Eqs. Q 
and ((BJ with the second Friedman equation ((2]), it is readily seen that the evolution of 
the scale factor for LJO model reads: 



2a-6. + h^ + k- SaHo'^a^ = 0, (7) 
which should be compared to: 

2afi + + A; - Aa^ = 0, (8) 

provided by the /ICDM model. The above equations imply that the LJO model has the 
same cosmic dynamics of a ylCDM Universe when we identify the creation parameter 
by the expression a = A/3Hq^ = f2j\. In particular, the LJO model predicts the 
same ylCDM transition from a decelerating to an accelerating regime with two basic 
advantages: (i) the cosmological constant problem is avoided, and (ii) the dark sector 
is reduced to a simple component (cold dark matter). The price to pay is the lack of a 
proper quantum field theoretical approach for irreversible creation of cold dark matter. 
Until the present, such a mechanism has been consistently justified only in terms of 
nonequilibrium relativistic thermodynamics 6, and kinetic theory S,. 

On the other hand, despite that both models (ylCDM and LJO) may share an 
identical Hubble expansion history, the same could not happen at higher orders on the 
theory of small density fluctuations. Therefore, it would be of interest to analyze the 
evolution of small density perturbations by taking into account the matter creation 
process. This is what de Roany and Pacheco proposed to analyze in their paper. By 
adopting the Poisson and Euler equations together a modified continuity equation (in 
order to include CDM creation) they performed a Newtonian analysis of the small 
density perturbations in the framework of the LJO model. However, as we shall see 
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below, the basic equations assumed in their analysis fail to recover the accelerating LJO 
model. In other words, even considering that the work is correct from a mathematical 
viewpoint, all the criticism in their paper refer to a decelerating cosmology and not for 
the scenario proposed by the authors. This is the basic result derived in the present 
note. At this point, we stress that the original notation of [2] will be adopted. In the 
notation of 1 , q is represented by Qv 

The Newtonian analysifl of de Roany and Pacheco [T] starts with a modified conti- 
nuity equation (in order to include creation of CDM) combined with Euler and Poisson 
equations: 

1^ + V ■ (pU) = pr, (9) 

^ + (U . V)U = -V<^, (10) 
^"^(j) = 47rG/9(t), (11) 

where p(t) is the CDM density, U is the velocity field of the fluid and F is the creation 
rate as defined in the LJO model (see Eq. ©). 

Let us try to see if the basic relativistic LJO equations are recovered from the 
above set of equations. To begin with, we recall that isotropy and spatial homogeneity 
of the unperturbed model implies that U = Hv (Hubble's law). In addition, since 
V ■ U = V • //r = 3//, one may rewrite Eq. (O as; 

p + 3Hp^pr (12) 

which is exactly equation Q of LJO. In addition, by inserting the expression of F as 
given by the above equation can directly be integrated. Apart a slightly different 
notation, the solution given by de Roany and Pacheco reads (see Eq. (13) in [1]) 

P ^ [Po - apcQ)a^^ + apco, (13) 

which also coincides with the solution of LJO (see Eq. (9) in [5]). Let us now consider 
the Euler equation (|10|l which can be rewritten as: 

Hr + Hr^{Hr)^{H + H'^)r = ~V(l>. (14) 

where the right hand side is defined by the Poisson equation (|lip . As one may check, 
since the density depends only on the time, an integration of the Poisson equation (|11|1 
yields 

V<^ = 47rGp^, (15) 

where the integration constant was chosen in such a way that the gravitational field is 
null in the center of the distribution. Finally, by inserting the identity H = a/a ~ 
in (|14|l and using (|15|l we derive the acceleration formula: 

a 4-kG , 
- = -— P, (16) 



^ It should be recalled that a nonrelativistic approach works when the scale of the pertur- 
bations is much less than the Hubble radius and the velocity of peculiar motions are small in 
comparison to the Hubble flow [9]. 
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which is fully different from Eq. ((5| of LJO model. This result means that the basic 
equations adopted by de Roany and Pacheco describe a decelerating model. In particu- 
lar, their description does not permit a transition from a decelerating to an accelerating 
Universe as required by SNe la observations. Naturally, for all practical purposes the 
note finished here. However, for completeness, it is interesting to know which is the 
Friedmann equation for the energy density in the kind of model discussed by de Roany 
and Pacheco. 

By inserting Eq. (|13|l into p6|l and multiplying the result by d, a simple integration 
yields 

.2 . , SttG 



a +k- 3 



(17) 



where k is an arbitrary integration constant. Finally, dividing both sides by and 
using again Eq. (|13p together ^ and it follows that 

„ k 



,=•2 

8ttG 



= 3^+3^, (18) 



which should be compared with the Friedman equation ([T]). In the above equation there 
is a spurious extra contribution, 12ttGpc- Its presence confirms again that the basic 
cosmological equations describing the unperturbed LJO model cannot be recovered 
with basis on the adopted nonrelativistic approach. As recently discussed by Basilakos 
and Lima TO', a consistent quasi-Newtonian treatment with pressure also requires a 
modification of the Poisson equation (see also [lip. 

In conclusion, we have shown that the Newtonian analysis performed by de Roany 
and Pacheco 1 fails when trying to recover the background (zero order) relativistic 
model proposed in Ref. [2]. All their results describing linear perturbations seems to be 
mathematically correct, however, they refer to a decelerating cold dark matter Universe 
with a modified continuity equation. It is also clear that by using different background 
equations will result on a different evolution for the linear density contrast since the 
corrections will generate new terms on the first order perturbation equations. A more 
rigorous "Newtonian formulation" for the relativistic LJO models, the corresponding 
evolution of small density fluctuations and other physical consequences is being devel- 
oped and will be published elsewhere |12] . 
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